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ABSTRACT 


The  problem  of  the  diffraction  of  electromagnetic  waves  by  a 
circular  aperture  in  an  infinitely  conducting  planar  screen  is  studied. 
Attention  is  focused  chiefly  on  the  low  frequency  end  of  the  spectrum 
where  the  product  ka  of  the  wave  number  k  and  the  aperture  radius  a  is 
small.   Transverse  electric  and  transverse  magnetic  modal  excitation 
appropriate  to  the  cylindrical  geometry  of  the  problem  is  assumed.   New 
integral  representations  of  the  solution  corresponding  to  incident  fields 
of  this  kind  are  presented.   These  representations,  which  involve  a  pair 
of  unknown  functions,  are  designed  so  as  to  automatically  satisfy  Maxwell's 
equations,  the  boundary  conditions,  the  radiation  condition  and  all  but 
one  of  the  edge  conditions.   The  unknown  functions  are  shown  to  satisfy 
a  pair  of  integro-differential  equations  which  follow  directly  from  the 
representations  by  means  of  an  elementary  calculation.   Convergent  power 
series  solutions  of  these  integro-differential  equations  are  readily  ob- 
tained for  small  enough  ka,  when  it  is  assianed  merely  that  the  excitation 
is  sufficiently  regular  in  the  aperture.   These  results  are  employed  to 
derive  the  first  tenns  in  the  power  series  expansions  of  the  far  fields, 
the  edge  fields,  the  aperture  fields  and  of  the  transmission  coefficients. 
When  the  incident  transverse  electric  and  magnetic  modes  are  specialized 
to  those  which  occur  in  plane  wave  excitation  these  expressions  are  found 
to  reduce  to  the  well-known  plane  wave  formulas  obtained  by  other  means. 


-Ill- 


ction  of  the  incident 
e  wave 


A  viev  of  the  diffracting  screen.  ^  is  the  screen,  ^  is  f^^fj^^^ll^^ 
of  radius  'a'  centered  at  the  origin  of  coordinates.  ^  +  OL  comciaes 
with  the  (x,y) -plane.  The  symbols  r,  ^  and  z  are  cylindrical  coordi- 
nates of  an  arbitrary  point  of  space;  r  denotes  the  distance  from  the 
T-ax  s  (th^  Lis  of  syLetry);  ^  is  measured  from  the  positive  x-axis 
in  the  X  v) -plane.  The  angle  0  is  measured  from  the  positive  z-axis 
to  the  ra/thSt  connects  thf  origin  to  (r,^z).  The  angle  y  fixes  the 
direction  of  propagation  of  the  incident  plane  wave. 
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1.  Introduction 

¥e  are  concerned  with  the  problem  of  the  diffraction  of  electromagnetic 
waves  by  a  circular  aperture  in  a  perfectly  conducting  planar  screen  when 

the  product  ka  of  the  wave  number  k  and  the  aperture  radius  a  is  small. 

r  1* 
Solutions  of  this  problem  have  been  given  by  J.  Meixner  [IJ  and  J.  Meixner 

and  W.  Andre jewski  [2]  using  series  of  spheroidal  wave  functions,  by  C.  J. 

Bouwkamp  [3a, b]]  and  N.  Marcuvitz   [3b]  using  formal  power  series  expansions 

in  ka,  by  H.  Levine  and  J.  Schwlnger  [k]   with  the  aid  of  variational  methods, 

and  by  Y.  Nomura  and  S.  Katsura  [5]  and  N.  Chako  [6]  employing  series 

expansions  of  hypergeometric  polynomials.   A  careful  survey  of  the  literature 

on  both  vector  and  scalar  aperture  problems  will  be  found  in  the  report  [3bJ 

and  paper  [^a]   by  C.  J.  Bouwkamp. 

In  this  paper,  we  present  an  approach  to  the  vector  aperture  problem 

which  is  an  extension  of  the  work  of  J.  Bazer  and  H.  Hochstadt  [7]  on  the 

scalar  aperture  problem.   Our  method  is  most  closely  related  to  two  recent 

attacks  on  the  problem  due  to  K.  A.  L-or'e  [8j  and  W.  E.  Williams  [9b]  '-^. 

All  three  approaches  ass^mie  excitation  of  E-mode  or  H-mode  type  such  as  is 

derivable  by  the  usual  operations  from  solutions  of  the  reduced  wave  equation 

having  the  form  cos  [m(f)-(f)^"^'^3u(  r  ,z)  ,  m  =  0,1,...;  here  r,z  and  ^   denote  the 

cylindrical  coordinates  of  an  arbitrary  point  of  space  [See  Figure  iJ  and 

(m) 
^  is  an  arbitrary  phase  angle.   In  addition,  all  three  approaches  have 

in  common  the  reduction  of  the  vector  problem  to  a  pair  of  scalar  aperture 

problems  by  means  of  the  Hertz  vector  potential  and  all  three  eventually  lead 

by  various  means  to  Fredholra  integral  equations  of  the  second  kind  whose 


^References  are  given  at  the  end  of  the  paper. 
**Footnotes  will  be  found  on  pages  kl   and  k2. 
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kernels  are  small  when  ka  is  small.   The  three  methods  differ  among  themselves 
chiefly  in  their  representations  of  the  solutions  of  the  scalar  problems 
mentioned  above  and  in  the  representations  of  the  vector  field  to  which 
these  give  rise.   The  essentially  new  element  of  this  paper  is  to  be  found 
in  our  representations  of  the  scalar  and  vector  fields,  which  possess  the 

following  noteworthy  features.   First,  our  representations  of  the  scalar 

h 


waves  have  the  relatively  elementary  form 


,  V  .  ,     /-,  --  \m  pi   ik  R  (r,z;t)   ,  r 

(1.1)    vi-"'  =  cOsCmH^']-"-  (i  ^     I         '-^^7—7-   4"''(*'"' 


^m  pi  ^ik  R  (r ,z;t) 

1    .,  „+ 


-^-1     R  (r,z;t) 


|m-l|  n      ^ik  R  (r,z;t] 

-1 
m  =  0,  1,  ... 


(1.2)    w^     =  cos[m<j)  -  ^^  ]rl    '    ;^       /   -^T " ^^   '(t)^t, 

R  (r,z;t) 


where  r"*"  =  v/r^+Cz+iat)^,  f^™^(t)  and  f^"^""'"^(t)  are  odd  and  even  functions, 

J--   1    jn      j_  .        n   n      J  T(ni)  iM  T(ni)   i(ni)   Jt 

respectively,  of  a  certain  general  class  and  9   =9    or     9^   =9   "^  o' 

Second,  when  ka  is  smalljthey  provide  the  basis  for  the  efficient  calculation 

of  explicit  formulas  for  the  field  everywhere  in  space  and  for  the  various  field 

characteristics  such  as  the  aperture  fields,  far  fields  and  transmission 

coefficients  for  all  m  >  0.   Third,  only  the  most  elementary  analysis  is 

required  in  dealing  with  these  representations;  no  use  is  made,  for  example, 

of  special  function  theory.   Finally,  the  field  representations  derived  from 

w-j   and  w^       have  certain  built  in  features  which  significantly  simplify 

the  task  of  meeting  the  conditions  of  the  problem.   Specifically,  the  field 

representations  are  designed  to  satisfy  automatically  Maxwell's  equations, 


the  boundary  and  the  radiation  condition  assuming  merely  that  the  unknown 

functions  flj   (t)  and  f^    (t)  exist.   Then,  the  condition  that  the  magnetic 

field  be  continuous  through  the  aperture  leads  directly,  by  means  of  an 

elementary  calculation,  to  a  pair  of  linear  integro-differential  equations 

for  f   (t)  and  f},    ("t)  of  a  type  already  discussed  by  Bazer  and  Hoch- 

stadt  [7]  in  their  treatment  of  the  scalar  aperture  problem.   Because  the 

Inhomogeneous  terms  of  these  integro-differential  equations  depend  on  an 

arbitrary  constant  (which  appears  linearly  in  both  terms)  the  reduction  of 

the  vector  problem  to  a  pair  of  scalar  problems  is  not,  however,  a  complete 

one  —  the  scalar  problems  are  coupled.   Nevertheless,  one  obtains  with 

relatively  little  labor,  convergent  power  series  expansions  in  ka  for 

f-j   (t)  and  f^  ~   (t)  when  ka  is  sufficiently  small  just  as  in  the  scalar 

case.   The  solution  of  the  problem  is  then  completed,  as  in  references 

[8] ,  [9]  and  [10] ,  by  adjusting  the  arbitrary  constant  so  as  to  ensure  the 

proper  edge  behavior.   It  may  be  mentioned  that  the  derivations  in  | 8 | , 

j_9a,  9bJ  ,  [10]  of  power  series  expansions  of  the  analogs  of  f-j   (t)  and 

f^  "   (t)  are  relatively  less  direct  than  the  one  which  we  have  just  described. 

By  the  method  just  outlined  we  have  calculated  for  each  m  >  0  the 
first  few  terms  of  the  power  series  for  f ^   and  f^     corresponding   to 
rather  general  classes  of  E-  and  H-mode  excitation.  We  have  nevertheless 
been  able  to  derive  for  each  m  >  0  the  first  terms  of  the  power  series  for 
the  electric  field  in  the  aperture,  the  edge  fields,  the  far  field  amplitudes 
and  the  transmission  coefficients.   When  the  E-  and  H-mode  fields  are 
specialized  to  those  which  appear  in  the  modal  expansions  of  plane  waves 


.1+- 


our  formulas  for  m  =  0,1  reduce  to  results  obtained  by  other  means  [see 
references  |j3a,  b]].  Our  formulas  for  large  m  seem  however  to  be  new 
although  it  should  be  possible  to  derive  them  from  Williams'  results  [9b] 
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2.   Formulation 

The  geometrical  setup  is  shown  in  Figure  1.  ^  and  -^  denote  the 
aperture  and  screen  respectively;   a  is  the  radius  of  the  aperture.   Together 
Wl  and  '^    make  up  the  x,y-plane.   The  origin  is  located  at  the  center  of  ^ 
and  the  z-axis  is  perpendicular  to  the  plane  of  the  screen.   For  the  most 
part  we  employ  the  cylindrical  coordinates  (r,(}),z)  , 

X  =  r  cos  (j), 
(2.1)  y  =  r  sin  cf),       0  <  (j)  <  2jt, 

z  =  z, 
r  denoting  the  distance  from  the  z-axis.   In  discussing  the  far  fields 
however,  it  is  more  convenient  to  employ  the  spherical  coordinates  (R,0,(})), 
where 

X  =  R  sin  0  cos  cj), 
,   ,  y  =  R  sin  Q   sin  ^,  0  <  (t"  <  2it, 

z  =  R  cos  0,  0  <  6  <  11^^ 


R  =  ]/r^   +   z^ 

The  angle  9  is  the  angle  between  the  positive  z-axis  and  the  position  vector 
to  the  point  (x,y,z). 

In  the  case  of  cylindrical  coordinates  e  ,  ei  ,  e  denote  unit  vectors 
at  (r,(f),z)  in  the  r,  ^   and  z  directions;  in  the  case  of  spherical  coordinates 
~R'  ~0'  ^l     denote  unit  vectors  at  (R,Q,([i)  along  the  radial,  latidudinal  and 
meridianal  directions. 

Assxjmlng  that  the  excitation  originates  in  the  region  z  <  0,  we 
represent  the  total  field  (  i^^ '^)  as  follows: 


z  <  0: 


(2.3) 


z  >  0: 


£(x)  =  E  (x)  -  E*(x*)  +  E*(x), 


^(x)    =  H  (x)  +  H*(x*)  -  H*(x*), 


Oi*    rw 


€(5)  =E(x), 
^(x)  =  H(x). 


Here,  (E^(x),  H^(x))  is  the  incident  field,  (E^(x)  -  E*(x*),  H^(x)  +H^(x*)) 

is  the  incident  and  reflected  field  that  would  be  present  in  the  absence  of  an 

aperture,  (E(x) ,  H(x) )  is  the  (unknown)  scattered  field  in  the  region  z  >  0 

and  (E-^(x*),  -H^(x*))  is  the  bacl-^scattered  field  in  the  region  z  <  0.   An 

asterisk  signifies  reflection  in  the  z-plane;  thus  if  £  =  (Q  ,  Q  ,  Q  ) 

denotes  any  Cartesian  vector  then  by  definition  Q-*  =  (O  ,  Q  ,  -Q  ) . 

*^  X       y  z 

We   assume  E  and  H-mode  excitation  of  the  form 


E-mode    (Transverse  Magnetic) 

(m). 


a)      E      ,    =  V       slnnr^        s  - 
~o,t  t  \  'in     oz 


met 


,(m) 


'm 


mcp   -   cp^    -  ,      (m  =  0,1,2,, 


(2.U)        b)      H^^,   =  ike^^V^    (sin4,,vf)). 


c)      E  =   -V?    (sinmA     V^""^)    , 

'        o,z  t    ^  "mo 


where    V^"^^    =V^"^^(r,z)    satisfies 
o  o 


1  a  /^^^o""^ 

r   Sr  \     Sr 


^2^(m) 


az^ 


0, 


k     - 


m 


V^"^)    =  0; 
o  ' 


-7- 


H-mode  (Transverse  Electric) 


^(m) 


a)   E   ,  =  ik  e  /\V,(cosm(i)  V^  -^ )  , 
~o,t      ~z    t^     T-m  o    ' 


(2.5)   b)   H  ^  =  -V^  (cos  mi    ^ 
~o,t     t  V     ^m   dz 


_  Sv^'") 


:(m) 


m<|)^  =  n4  -  <|)^"'^    (m  =  0,1,2,...) 


c)   E^  .  =  0, 


H    =  V^(cos  mi  V^), 
o,z    t      ^m  o  ' 


where  V 


(m)  _  -Jm) 


=   V   (r.z)  satisfies 


1  a 


av 


,v(«i)\    ;^2-( 


a"v 


,m 


^  V   ^ 


+ 


2   m 


V 


(m) 


In  these  equations  (j)^   and  ^         are  arbitrary  phase  angles,  V 
is  the  operator  V  -  e  -r-      (V  is  the  gradient  operator  in  three  di- 
mensions) and  Q,  =Q-eQ,Q=(Q,Q,Q)  being  any  vector. 
~t   ~    z  z  ~     iC  y  z 

We  represent  the  scattered  field  in  terms  of  the  Hertz  vector 
potential  II, 


(2.6) 


E  =  ik  V  /N.  n 


H  =  k  n  +  v(v.n), 

where  (in  Cartesian  coordinates)  II  satisfies 


(2.7) 


v^  n  +  k^  n  =  0. 


Assuming  E-mode   incidence   and  setting    [cf.    references    [8] , [9 | , [lO] ] 
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(2.8)  n^  =  -Sin[(m-1)(^-<|)^"^)]V^"^-I)(r,z), 

n^  =  cos[m<j)-<{)^"'^]v["'^  (r,z),  m  =  0,1,2,..., 

where  V^   (r,z)  and  V^  ~   (r,z)  denote  two  unknown  scalar  functions  of  r  and 
z  we  find  immediately  that  the  radial  and  (|)-coiiiponentes  of  II  are 

n^  =  cos  m(})^  V^"""  ''  (r,z),         mljj^  =  mcj)  -  ^^^' 
(2.9) 

n,  =  -sin  mc})^  V^~      (r,z),  m  =  0,1,2,... 

from  which  it  follows,  owing  to  equation  (2.6),  that  the  corresponding 
components  of  E  are 

E   =  ik     2 .  21  yW  )    sin  mi   , 

r      \    oz      r     1       f  ^m   ' 

To  obtain  the  corresponding  H-mode  formulas  we  have  only  to  replace  V^    , 

Y^^\   ^l"^^    and  (f,^  by  Y^^'^K   Y^^\  Pj"^    and  ^^  in  equations  (2.8)-(2.10) 

(see  footnote  7) . 

Suppose  now  that  we  are  successful  in  finding  functions  n  ,  IE  and  n 

X   y     z 

of  the  form  of  the  right  members  of  equation  (2.8)  having  the  following 

properties: 

(i)  n  ,  II  ,  n  satisfy  the  reduced  scalar  wave  equation. 
^     x'  y'  z 


p.v(m-l) 


(ii)      Screen  conditions 
Sv(m-l) 

(2.11)  0  =  -A —  =  vy 

oz       1 

(ill)   Radiation  conditions;   as  R  approaches  infinity 
I    \  (    \  ikR 

vf)~  a|-)(o)  s__, 

(2.12)  .^ 

v(m-l).A(m-l)(e)  e_.^  „,  =  0,1,2,, 

(iv)   Edge  conditions;   V-|   and  V^     are  such  that 
Uni-l)   3  (m) 


(2.13) 


in  the  aperture  {X  at  a  distance  p  from  the  circular  edge. 

(v)   The  aperture  condition.   If  the  excitation  is  of  E-mode  type 
then  in  the  aperture  (x  {^  <   a) 

3vjm)  ^        _       ^2^^m) 

(2.li|) 

„(m-l)  /,     vni-l        .,  -1   -m  S      /  m  ,r(m)x 

YX  =   -mc    (kr)  -  ik     r       -sr  (^     ^       ) 

2  m^  or    ^  o 
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If  the  excitation  is  of  H-mode  type  then  in  the  aperture  (^  (r  <  a) 

(2.15) 

(m-l)  -, 

V^    =  mc  (kr) 
2       m 

In  equation  (2.1^+)  and  (2,15)  c  and  c  are  arbitrary  constants  vhich 
will  later  be  employed  to  satisfy  the  edge  condition  (iv). 

We  maintain  that  if  such  functions  can  be  found,  then  we  have  in 
effect  solved  the  vector  problem.   For  if  condition  (i)  is  satisfied  then  H 
satisfies  the  equation  (2.7)  and  the  derived  quantities  (E,  H)  of  equation  (2.6) 
satisfy  Maxwell' s  equations.   As  a  mere  inspection  of  equation  (2.10) 
reveals,  (li)  and  (iv)  are  simply  conditions  which  are  sufficient  to  ensure 

Q 

that  E^  vanishes  on  -^  and  has  the  proper  edge  behavior  .   Furthermore,  the 
~t 

radiation  condition  (iii)  and  equation  (2.6)  imply  after  a  brief  calculation 
that  as  R  ->  00 


(2-16) 


where 


(2.17)  fS^ho-.^J   =   (A^""-^'    sin  0  +a[°'   cos  0)    cos  m^^  g^ 


+   (a'""-"-'    cos   0   -  a'""'    sin  0)    cos  mcl     e 


m  ~© 


A (m-l)    .        i 
-  AA  sm  md)     e  i 

2  ^m  ~<p 
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lt follows   that  R|E|    and  R|H|    are  boimded  when  R  ^  <»  and  that 
(2.18)        0  =     lim     [r[^^^E-H]J     =     liiii     (r[^j^-^  H   -  Ej] 

R->ooL,  J  R->ool^  J 

which  are  the  appropriate  radiation  conditions  for  the  electromagnetic  field 

[see  reference  Ll^JJ.   Finally,  as  we  show  in  Appendix  I,  the  aperture 

conditions  are  equivalent  to  the  requirement  that  7/  be  continuous  through 

the  aperture.   It  is  a  well-known  fact  that  this  condition  guarantees  the 

continuity  of  Q^   and  the  first  derivatives  of  (c,  n^    through  the  aperture 

provided  that  //  is  sufficiently  differentiable  in  Ot. 
~t 

For  technical  reasons  of  this  kind  it  is  necessary  to  place  restrictions 
on  the  incident  field  in  the  aperture.   The  following  conditions  are  sufficient 
to  ensure  the  success  of  the  ensuing  analysis  but  they  are  by  no  means 
necessary  conditions.   Assuming  E-mode  incidence,  we  require  that  V   , 
SV   /Sz  and  S  V   /Sz  have  the  following  form  in  the  aperture  Cc. 


H  =  0   ^ 


vf^  =  (kr^  YL       h^°^(m)(kr)2^ 


(2.19)  ^  =  k(kr)"'  ^   h^^^(m)(kr)^^,         (E-mode) 


H  =  0 


32^(m) 


^-  =  ^{^tY-     h^2)^)(j^)2^^ 


az2  ^  =  0  ^ 


and  that 


9 


(2.20)  ^2_  =  l,p(m)^(m)  ^.^^^ 
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where  p    is  a  constant.   The  three  series  are  assumed  to  converge  in  a 
circular  domain  of  radius  a  -f  A^  A  >  0  centered  at  r  =  0;   the  coefficients 
h^      ,  j   =  oP-,2;  M,  =  1,2,..,  are  assumed  to  be  independent  of  k  but  may 
depend  on  m.   The  corresponding  requirements  for  H-mode  incidence  may  be 
obtained  from  the  above  simply  by  replacing  V   ,  h^"^   and  3^   by  V^   , 
h^"^   and  3^   ,  respectively. 

Suppose  the  incident  field  is  a  plane  wave  whose  direction  of  pro- 
pagation n  is  in  the  (x,z) -plane  and  makes  an  angle  7  <  rt/2  with  the  positive 
z-axis.   Such  an  incident  field  may  be  expressed  as  a  superposition  of  the 
plane  waves 


(a; 

(2.21) 


E  =  -e  ^  n  k  exp(ik(n-x) ) ,     H  =  -e  k  exp(ikn-x): 
~o    ~y   ~         ^~  ~' '  ^     ~o    ~y         ~  ~ '  ' 


(b)    E  =  -e  k  exp(ikn-x)  ,  H^  =  +e  /^  n  k  exp(lkn'x) 

where  e  is  the  unit  vector  parallel  to  the  y-axis.  Using  the  Fourier -Bessel 

~y 

series  expansion  of  the  first  factor  of  exp(ikr  sin  7  cos  (|))exp(ikz  cos  7")  = 

=  exp(ik-x)  we  find  that 

•  m  ^   -P  /-,    .    N   ikz  cos  7 
/  N     ,    s      1-     K     J   (kr  sm  7)  e 

mj    T^i"^)  m  m ^ 


(2.22)    V^'''^  =  V 


o      o  sm  7 

(k     =  1   when  m  =  0:   K  =  2  when  m  =  1,2,. . . ) 
m  m 

and  that  for  all  m,  including  m  =  0,the  products   (mc})  )  and  (m^   )  have  the 

values 

(2.23)         (m^J    =m^+  f,  i^J    =   4- 
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The  functions  V^"^  and  V    evidently  satisfy  all  the  requirements  listed 
o        o 

above  and  from  the  power  series  representation  of  J  (kr  sin  7)  we  find 

easily  that 

1   ^   K      (sm  7) 


,(0),  \   r(o),  >,  m 

h  (m)  =  h^   (m) 


2     i-il  r(m-tiJ.+l} 


i^\  .m+2|J.+l  ,^  ,  .    an+2^-1 

(1)  _, -.X  X       ^        K     COS  7  (sm  'O 

(2.25)   h   (m)  =h^"^m)  =          "^ 

^  ^  'f'^^^   n:  r(m-Hx+l) 


.m+2iJ.+2  ,,     2   ,  .    ^m+2n-l 
1       /<  cos  7  (sm  7j 


A2),    .        r(2)/  \      m 

h    (m)  =  h    (m)  = — — 

^       ^  ^^^^   n;  r(m-Hx+l) 


and  that 


(2.25)         p^'"''  =  p^"'''  =  i  cos  7 
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3.   Representations  of  v["^^  and  y^"^"-*-^  -  E-Mode  Excitation?-'"' 
3.1  The  Representations.   In  this  subsection,  we  describe  the 


representations  of  V-[   and  V^  ~   associated  with  E-mode  incidence  and  prove 
in  detail  that  conditions  (i)-(v)  of  Section  2  can  be  satisfied.   The 
corresponding  analysis  for  the  case  of  H-mode  incidence  is  completely  analogous 
and  will  therefore  be  omitted. 

Following  J.  Bazer  and  A.  Brown  [lla,bj,  W.  D.  Collins  [12]  and 
J.  Bazer  and  H.  Hochstadt  [?]  ,  we  represent  Y^'    and  V^  ~   by  the  expressions 

/  N     I  I  /  1  :n  \|ml  p    ikE  (r,z:t)   ,  s 

(3.1)   v"^)  =-'"'^|:    r  ^  /     ff^(t)dt 
J         V^  *y    J__^         R  (r,z;t)    ^ 


1    ^    .,^+ 


m 
r '  ' 


i  -2-      /  <  ^— +  (-1)^  ^ >  f^."'^(t)dt 


"■  ^^     i     \      R^Cr^zjt)     "      R-(r,z;t)   '   ^ 


for  j  =  1,2.   where 


(3.2)         lf(r,z;t)  =    /r^+(z  +  iat)^ 


Hereafter,   it  is  to  be  understood  that  m  may  be  any  non-negative 
integer  when  J  =  1  and  that  m  =  -1  or  any  non-negative  integer  when  j  =  2. 
To  fullfill  the  requirements  on  V-[   and  V^   we  assume  that  the  f  .   (t) 
have  the  following  properties: 

(i)   f .   (t)  is  even  or  odd  according  as  j  is  even  or  odd;  we  have 

J 

already  made  use  of  this  property  in  the  second  expression  for  the  V.   . 

J 


-15- 


(ii)   f-,   (t)   and  f^   (t)  are  regular  in  t  in  the  circular  domain 

defined  by  the  relation  |t[  <  \+i\/a,   A  >  0  for  all  sufficiently  small  values 

12 
of  ka  and  all   m>  1.  Since  the  regularity  of  the  incident  field  in  the 

domain  a  +  A^  A  >  0  [see  the  remarks  following  equation  (2.20)3  ultimately 

imply  that  f .   (t),   j  =  1,2  have  this  property,  no  generality  is  lost 

J 

in  assuming  it  at  the  outset.   Finally,  to  ensure  the  proper  edge  behavior 
[see  equation  (2.13)]  it  turns  out  to  be  sufficient  to  require 
(iii)  that  for  each  m  >  1   (D  =  t— ) 

,    ^     dM"^^(1)  =0,  k  =  0,1,...,  m-l' 

(3.3)      ^ 

pk^(m-l)^^^  =  0,  k  =0,1,...,  m-2 

and  that 

(3.^)  dI"^-1|  f^'^-l^d)  -  ^-  D^  f^)(l)  =  0. 

When  m  =  0  condition  (3-3)  is  to  be  ignored. 

To  resolve  the  ambiguity  in  the  definition  of  R  (r,z;t)  when  z  =  0, 
we  assume  that 


(3-5)         lim  R-(r,z;t)  =  \/r^  -  a^t^  ,         when  r^  >  a^t^ 
z  I  0 


+     /~2     _    2 
which  implies,  on  choosing  the  branch  of  R  =  \/r  -(at  +  iz)   in  an  appropriate 

manner ,  that 

(3-5')         lim  R-(r,z;t)  =  ±l\Ja^t^   -   r^  ,       when  r^  <  a^t^. 
z  i  0 

At  points  in  the  aperture,  (z  =  0,  r  <  a)  the  only  points  in  physical 

space  where  the  integrands  of  equation  (3.1)  may  be  singular,  V.   ,   j  =  1,2 

J 

and  its  first  and  higher  order  derivatives  are  defined  to  be  the  limits 


-16- 


approached  when  z  approaches  zero  through  positive  values.   The  limits 

themselves  can  be  rigorously  derived  with  the  aid  of  the  properties  of  the 

f    \ 
f .   (t)  listed  above  and  standard  function-theoretic  techniques.   Appendix  I 

J 

of  reference  [llaj  relating  to  the  calculation  of  the  edge  behavior  exemplifies 

13 
in  detail  the  sort  of  argument  required.    Here,  it  is  enough  to  say  that 

the  several  limiting  processes  employed  in  what  follows  are  justifiable  and 

that  their  end  results  agree  with  those  that  are  obtained  by  purely  formal 

considerations . 

If  it  is  assumed  merely  that  the  f':™^(t),  j  =  1,2,   satisfying  (i)-(lii) 

J 

above  exist,  then  it  can  be  shown  even  at  this  stage  that  v|   and  V^ 
satisfy  all  the  requirements  listed  in  Section  2  with  the  exception  of  the 
aperture  condition.   We  must  show  first  that  11  =  cos  mcj)  V^'^'^  (r ,z)  and 
(n^,  n  )  =  (cos(m-l)<|)^,  -sin(m-l)^^)V^"'"^^r,z)  with  Y^^^    and  V^"""^^  defined 
as  in  equation  (3.1)  satisfy  the  scalar  wave  equation 

(3.7)  2^  .  iff  .  ^  .  k==n  =  0. 

^    by         Sz 

The  proof  is  essentially  the  same  for  11  ,  II  and  n  so  that  it  will  be 
■f  >'  x   y      z 

sufficient  to  explain  it  in  relation  to  one  of  these  components  —  H^  say. 
Wow  let  W(x,y,z)     be  any  solution  of  equation  (3-7) •   Then  cW/Sx  and 
Sw/Sy  are  also  solutions  and  hence  so  are  the  functions. 

(3.8)  .^  =  (|-.i|)„. 

In  fact, 

(3.9)  (D*)n,  =  (I;  ±  f  f  W 
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are  solutions  of  equation  (3-7),  m  =  1,2, Expressing  D  in  polar  co- 
ordinates [see  equation  (2.1)^  we  arrive  at  the  operator  identity 

(3.10)         B*  =  e±^+(^  i  i|)  . 

If  D~  is  now  applied  to  any  axially  symmetric  solution  W  =  U(r,z)  of  the  wave 
equation  it  follows  easily  that 


(D^)-U 


tun 


4 


S_   m-l 
Sr    r 


B_   m-2 
Sr    r 


a_  _  1 

(^   r 


Br 


U(r,z) 


which  may  be  written  more  compactly  as 

(3.11)         (D=^A(r,z)  =e±^-^r-(i|:ru(r,.z). 

Evidently,  any  linear  combination  of  (D  )  iJ  and  (D  )nj  is  a  solution  of  the 
scalar  wave  equation  provided  U  is.   If  we  set 

P        ikR'^(r,z;t) 


(3.12) 


U(r,z)  = 


(^   R^(r,z;t) 


f^.^)(t)dt 


then,  since  U  is  obviously  a  solution  of  the  scalar  wave  equation,  so  is 
the  following  linear  combination  of  (D  U)   and  (D  U) 


(3.13)     n  =  cos 


/  \.   /  1  :i  \™  P    i^  (^^z;t) 
L  H^   ^  J   Vr  Sr/  j    R+(r,z;t) 


f^."^^t)dt. 
J 


Condition  (i)  of  Section  2  is  thus  satisfied. 

To  verify  that  h"^]^'^' /hz   and  vj"^^  vanish  on  <^ [condition  (ii)J 
we  need  only  note  that  the  integrands  in  the  last  member  of  equations  (3.1) 
and  their  first  derivatives  are  non-singular  when  r  >  a;  the  vanishing  of 
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SVp  "  /5z  and  V^   then  follows  immediately  on  letting  z  approach  zero 
under  the  integral  signs. 

To  verify  the  radiation  condition  (condition  iii)  ve  make  use  of  the 
spherical  polar  coordinates  of  equation  (2.2)  and  observe  that 

R"(r,z;t)  =  R  ±  lat  cos  9  +  Q'i'R''^) 

when  R  is  large.   With  the  aid  of  this  formula  it  follows  immediately  that 
equations  (2.12)  are  satisfied,  a'"^^(0)  and  A^  ~   (6)  being  defined  by 

1 
(3.1U)         A^^\q)    =   -2(ik  sin  0)""  rsinh(atcos  0)f [""^  (t)dt, 

0 


(3.15)        A^'^'^^Q)  =  2(ik  sin  0) '"""^1  T  cosh(atcos  0)f2^"'"^^  (t)dt 


To  prove  that  the  edge  conditions  (iv)  of  Section  2  are  satisfied, 
it  is  necessary  to  determine  the  behavior  of  the  functions  Sv),    /Sz, 
V     and  Sv_[   /^  in  the  neighborhood  of  the  circular  edge.   To  exemplify 
our  procedure,  we  show  how  the  behavior  of  V^    is  obtained.   We  begin  by 
assuming  that  | z/a |  <  e,  where  e  <  r/a  is  a  small  positive  number.   Then  we 
deform  the  contour  near  t  =  0  in  the  right  member  of  equation  (3-1)  into  a 
semi-circle  of  radius  e  centered  at  t  =  0  and  falling  below  the  real  t-axis. 
On  the  new  contour  --  C  say  —  the  ratio  | z/a |  will  then  be  everywhere  less 
than  |t|.   Then,  we  observe  that  the  integral  representation  of  V.^    [see 


equation  (3.1)^  niay  ^^   expressed  as 


(3.16) 


V. 


(m) 
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ir\      j       expiik  R   (r  ,z  ;t) ] 

^^  R   (r,z;t) 

e 


1 


dt  Vz+iat 


m 


4"^(t)dt, 


on  noting  that 


(3-17) 


1   S_   fexpji'kR'^) 
r   Sr    V       p+ 


ia(z+iat)    dt 


d       / eyp(ikR 


R 


+ 


integrating  by  parts  m  times  and  finally  making  use  of  the  "behavior  of 
f|"^''(t)  at  t  =  ±1  [see  equation  (3. 3),  and  recall  that  f|™^(t)  is  oddj  . 
Next,  we  make  use  of  the  fact  that  exp(ikR  )/r  =  (z+iat)   d(exp  ikE  )/dt, 
integrate  by  parts  once  again  and  find  that 


,+/ 


f3  18)    V^™^  =  -  —  l—]'^   rexpfikR""(r,z;t)] 
^^'      '  1       ka  Va  /  I     z+iat 


dt  \ z+iat/ 


m 


f^"^)(t) 


-1 


expfikE  (r,z;t)] 
R"^(r,z;t) 


dt  V  z+iat 


-im+l 


f{"^)(t)dt. 


If  we  set  z  =  0  and  write 


(3.19)        r  =  a-p, 

.+   .    /.,  ^+\  . 


p  «  a 


1^1  . 


and  expand  R  and  exp(ik  R  )  in  powers  of  p/a  we  find  that    in  the  aperture  (^ 

.  d"^f^)(l)     r 


(3.20) 


.(m) 


1  .    ra+1  ^m 

la  dt 


"V^(p/a)l/2   +Cr[(p/a)3/2j 
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Similarly,  we  find  that        in    /^ 
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(3.21) 


(m) 


^(m)^(in)^^)  r 


m+2 


la 


dt 


m 


x/2(p/a)-l/2  +  (7[(p/a)^/2j 


(3.22) 


dV, 


(m-1) 


1/1  X  ^^-'^-'(^ 


dz 


a  \  m 
-a 


df 


m-1 


-\/2(p/a)'l/2  ^  ^^^p/^)l/2j 


It  can  also  be   shown,  in  the  same  way,  that  [see  reference  [tIJ 


(3.23) 


V 


m  (m) 
(m)  _  2V2    ''^l)   ,./.^l/2  ^, 


1      .   m+1     ,,m 
1  a       dt 


(p/a)^/^  sin  (6/2)  +  (7  [(p/a)] 


(3.21^) 


r(m-l)  ^  y(hi-1 


-2\/2 


:r  /I  \  d 


r=l 
z|0 


a  /  dt 


-j^f?-'(l)(p/a)^/2^os|-^S[(p. 


Here,  we  have  put 


(3.25) 


z  =  p  sin  8, 

r  =  a  +  p  cos  5, 


[see  Figure  2].   V^"^""^^ 


0  <  S  <  rt. 


r=:l 
zjO 


is  a  constant  whose  value  may  be  calculated 


from  equation  (2.l4)  once  the  constant  c  has  been  fixed.   V.;   and  y\ 
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in  equations  (3.2^1)  and  (3.25)  may  be  differentiated  formally  with  respect 
to  p  and  6;   the  resulting  formulas  then  lead  to  the  correct  expressions 
for  (E,H)  in  the  neighborhood  of  the  edge. 

Before  outlining  how  f-j"^  (t)  and  fp  ~   (t)  are  determined,  we  should 


mention  that  there  are  several  ways  to  represent  the  solution  of  the  scalar 
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problems  formulated  in  Section  2.   In  addition  to  the  works  already  cited 
at  the  beginning  of  this  section  mention  should  be  made  of  the  works  of 
D.  S.  Jones  [l3a] ,  A.  E.  Heins  and  R.  C.  MacCamy  [15]  and  ¥.  E.  Williams 
[pbj  .   How  these  methods  are  related  to  each  other  is  discussed  by  B.  Noble 
in  [16]    and  by  J.  Bazer  and  H.  Hochstadt  in  the  appendix  of  [7 1 . 

3.2  The  Aperture  Conditions  and  the  Integral  Equations. 

We  begin  by  computing  the  limiting  value  of  Sv^  /Sz  in  the  aperture. 
Assuming  z  >  0,  we  differentiate  the  first  and  last  member  of  equation  (3-1) 
and  find  for  j  =  1  that 

(3.26)        ^  =  -"  (f  hT"''  f    ^^^^^^^    e^l±^'-{r,.-t)]   ^(m)  ^^^^^ 

0  ^  (r.z;t) 


m  fl  b  \  r    ,      .  ^v  expfikE  (r,z:t)]  ^(m),J_^-,^ 

-^  [r^l        J      (2-iat)  --t^— — ^^--^  f|  ^t)dt; 


owing  to  the  identity 

f^  071  ^    /expfikR-]\ 1        1   a     /expfikR-] 

^^      '^  ^V         ^±  y~z±iatr^V        o± 

n      '  n 

Letting  z  approach  zero  through  positive  values,  and  making  use  of  the 
limiting  behavior  of  R"(r,z;t)  ["see  equations  (3.5)  and  (3-5' )J  we  find 
that  in  W 
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^  =  2i     i|-  r        — === tfW(t 


(3-28)  -""-^=2i^-^j         ^  r-     2      2 

0  \/(r/a)^-t^ 


)dt 


,m+l    p''/''   sin(a\/(r/a)2-t2) 


-2(i^)       [      -     , — -        tfM(t)dt 

;r/a)==-t' 


^   Set/         j  ^-22 


1  /  p              p 

,^    ;.    vin+l     n  sinh(aV't    "(I'/a)    )           /    ) 

-2-1:  ^== *  f      ^^)<^* 

V^    ^/         J,  ,  A2    ,    /    ^2                        ^ 

r/a  \/t   -(r/a) 


Since   sin(a\/(r/a)2-t2)    /\/{^af-t^    ,   0  <  t^  <   (r/a)^,    is  the  continuation 
;aV^^-(r/a)V\A^-(^/a)^     when  t^  >   (r/a)^,  we   conclude  that  for 


of       sinh 
0  <  r  <  a 


/S'   oo<=('„NAWn^2    ,2 


m+1    p'      cos{(x\/{T/a)    -t      ) 


,    :s  \m+±    p'       cosi^av  V-c/a';    -'-      -^  Cml 

(3.29)      li|         /       -=== *4"'(*'« 

0  V{r/a)=^-t2 


1  /        p  P 

1  /l   S  \"^'*'"^    r      sinh  (a  V*   "(^/a) 
0  Vt^-(r/a) 


i  i  4-       /  =— -— = tf;"'^(t)dt 

i   Vr   dry         J  r-^ — ^  1 


,-m     dvj^'^r^O) 


21  ^^ 


Wow  from  equations    (2.lU)    and   (2.19)    it  follows  that   Bv|"^'' (r,0)/5z, 
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rW 


(3.30) 


M.=0 


where  h  (m)  is  independent  of  k.   Introducing 


(3.31) 
(3-32) 


T]  =  (r/a)   , 

e  =  t2 


into  equation  (3.29)  and  then  setting  ^  =  t]  | '  in  the  first  integrand,  we 


find  after  dropping  the  prime  on  I '  that 


(3.33) 


dTi) 


m+1 


(°^)  f  ,  /7^ 


n/^  ff^  (\/ti?  )cos[a\A(l-0]   ff^\/l)sinh^\/Ti-q 


Vi^ 


i\/Ti-l 


dl 


m+1  m+1  v —        ^ 
a   a  )  ^,(■^/''^^l-'• 
^^:;5—  ^  h  (m)(aTi) 

i2  |J.=0   ^ 


Proceeding  in  an  anologous  fashion  with  Vp   (r,z)  we  note  from  equations 
(2.IU)  and  (2.19)  that 

r(m-l) 


(3.3^)   V) 


,2n 


,  m-1  |m-l|  \     /  w,  N^M- 
=  k   r'    '  /_  p  (ra)(kr)  "^ , 

z|0  n=0  ^ 


0  <  r  <  a 


where  p  (m)  is  independent  of  k  and  find  that  f^  ~   satisfies  the  following 
integral  relation 


(3.35) 


m-1 


dri 


,(m-l),.  /-rr^-.r...  /TTT-m        f(m-l), 


f  2  "  ^ '  (  \AilTcos [a\/Tl(l-Oj        4   ^^(\/g)sin]ri(a0H: 


V^  Vi-^ 


i\/i    \/r\-i 


d? 


m  m-1 


a 


■  ™~       \ —  2 


M.=0 
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for   m  =  1,2, . . .  . 

Equations  (3.33)  and  (3-35)  s.ve   the  integral  relations  mentioned 
in  the  introductory  section.   Note  that  they  follow  from  the  fundamental 
representations  (3-1)  in  essentially  one  elementary  step. 

Since  equations  (3-33)  and  (3.35)  are  basically  of  the  same  form  it 
vill  suffice  to  explain  how  equation  (3-33)  is  solved  for  f^      (t). 

¥e  begin  by  assuming  that  a  is  small  and  seek  a  solution  f.:   (t)  of 

equation  (3.33)  having  the  form 

00 
(3.36)        f^^^^t)    =  t(l-t2)-  a"^-^^  ZI     B^-)(t)a^ 


1/2/-,      X   m+1  y~        {m),r-.    \x        „(m)/    /-. 
=  T)         (l-Ti)a  /__     B'      (Vil)cL     =  f-,       (\/ti) 

[1=0        ^ 
(t]   =  t     here) 

where  B   (t)  is  a  polynomial  which  is  even  in  t  of  degree  at  most   |a. 

It  has  been  demonstrated  by  J.  Bazer  and  H.  Hochstadt  |_7,pp.  17-19  |  that  a 
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solution  of  this  form  exists  and  is  unique  for  all  sufficiently  small   a, 

the  power  series  converging  for  all  t  such  that  |t|  <  1  +  A/a  [see  the 

remarks  following  equation  (2.20)  and  preceding  equation  (3.3)].   This  being 

the  case,  we  may  substitute  the  expression  for  f^      (  \/t )  into  equation  (3-33) j 

expand  in  powers  of  a  and  equate  coefficients  of  like  powers.   We  then  find 

that  we  are  able  to  calculate  the  B^  -^(t)  recursively.   To  convey  the  basic 

idea  of  the  procedure  let  us  find  B   .   Since  the  degree  of  B    is  at 

o  ^        o 

most  zero,  BW  must  be  a  constant.   Equating  coefficients  of  a    we  obtain 
from  equation  (3-33) 
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(3.37: 


(m)  d^::^  r  n(i-^ir  ^^  ^ 


m+1 


o    ,  m+1 
dt) 


yrn 


i  2 


m+1  o 


h  (m) 


owing  to  the  fact  that  the  term  involving  sinh[ayTi-|  ]  is  of  order  a 
From  this  relation  it  follows  directly  that 


m+2 


(3.38)  B 


,    ^      /    -,  sm     m+1  h   (m) 
(m)_   (-1)      a o^ 


°     ~  2^'-^i^n     r(m-^|) 


The   corresponding  expressions  for  B^      (t),   B^      (t),   B,       (t)    are 


B^^t) 

- 

0, 

B^"^)(t) 

= 

(-l)"^a'^^^ 
iV^2"^+^(m+|) 

h^(m) 

hjm) 
Mm+  |)  _ 

B^"^)(t) 

= 

2ah^(0) 
-       „  2        \o' 

+  t 


h^(m)+  -j^ 


(3.39) 


m 


(m+l)h   (m)  2m2_^+, 

(m+l)h   (m)- T +      ^     '?^^ 


■A 


i  y^2"^-^^(m4l) 


\  +  mt 


2  Mm  +  i)  6i|(m-|)(m+|)2      ° 


2h2(m)    +±h^(m)    -  — 

I6(m  +  :j) 


h^(ra) 


V 


t^  2h2(m)    +  i  h^(m)    +  ^  h^(m) 
Similarly,  we  assiome  that  f^  (t)   has  the  form 


(3.^0) 


A'^-^ht)    =    (l-t2)l"^-l|    a"^-l  21  C^"^-^^t)a^ 
■^  M.=0      ^ 
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where  C^"^   '^(t)  is  an  even  polynomial  of  degree  at  most  |-i  and  find  that 


equation  (3-35)  implies  for  m  >  1 


I    -,  ^m-l  m   ,  ^ 


m-1 


\A  2--  r  ( 


m 


1, 


^"'"^^t)  = 


2ap^(l) 
_ 

rt  1 


Ira 


4"-^\t: 


/  , ^m-l  m 
(-1)   a 


r 


sA2^-lr(m  +  |) 


(m-1) 


P-L(m) 


Pjm) 


Mm-|) 


+  f 


p^(m)  +  ]j;PQ(m) 


(3A1) 


^a  p^(l) 

3    ^M  ' 


c^'^-^^t; 


2a 
2. 


'^o(i)(|-^)  +  fpid)  +  ^  Po'i) 


Im 


^a^p  (2) 

^~2^    ^2m  ^ 
9rt  1 


and  for  m  =  0  that 


(3.te) 


H  =  0,1,2,. 


It  is  to  be  understood  here  that  the  p  (l)  which  appear  in  the  expressions 
for  the  C^   (t)  in  equation  (3- ^2)  are  to  be  replaced  by  p  (O).) 
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3.3  Ensuring  the  Proper  Edge  Behavior 

As  a  mere  glance  at  the  form  of  f   (t)  and  f^  ("t)  reveals 

|_cf.  equations  (3-36)  and  {3.hO)J ,   the  edge  conditions  (3-3)  are  satisfied. 

¥e  have  yet,   however,  to  satisfy  condition  (3.''+).   To  this  end,  ^t . 

references  [8]  ,  [9I  j  D-'-jJ  ^^  recall  that  in  the  aperture  C(  [see   equation 

(2.14)J 

^M  ^2^(m) 

1  .  c  k(krr   ^    ° 


dz     m  '   '    k   ^  2 

dz 


(3.^3) 


,r(m-l)       /,  sxa-l        .,  -1  -m  S  /  m  „(m) 
VA     =  -mc  (Icr)    -  ik   r   -^r-  i"  V^ 
2         m  or  V    o 


and  note  that  each  term  in  the  right-hand  members  of  these  equations  is  of 

the  form  assumed  in  equations  (3.30)  and  (3.3^)j  respectively,  and  that 

c  is  an  arbitrary  constant.   Next,  we  write 
m 

{3.hk)  f["^)(t)  =c^4^{  (t)  +fj^)(t) 

(3A5)       f^^-'-ht)   =   c^  f^^-^^t)  +  f^^-'^^*)  ' 

where  we  suppose  that  f\^\    (t)  and  t\^l{t)    are  solutions  of  (3.29)  with 
k   /2i  and  (-S  V   /Sz  )/(2k  r  )  as  inhomogeneous  terms;  we  define 
f^  ^   (t)  and  f^  ~  '^(t)  with  respect  to  the  analog  of  equation   (3.29) 
in  a  similar  manner.   In  this  way,  we  obtain  four  functions  f .[  . ,  f ^  .   , 
j  =  1,2  with  the  following  properties:   i)   f|'^|(t)  ,  f|  gft)   have  the 

same  form  as  the  right  member  of  equation  (3. 36)  while  A^~    ^(t)  and  f^"^""^^  (t) 

2,1  '^}'- 
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have  the  same  form  as  the  right  member  of  (3'^0);   ii)   f^n  -("t)  and 
(m-1) 


f^"\'^'(t),  j  =  1,2,  are  independent  of  c^. 


IS 


22 


To  satisfy  the  edge  condition  (3.k)   we  now  substitute  the  above 

Jm-1), 
"2 


expressions  for  f;"^^(t)  and  f'"'  ^\t)    into  equation  (3.^)  and  solve  for  c  .  The  resul 


"1 


m 


(3.U6) 


c  = 

m 


-iaaD'"-^'f^";-^^(l)  •^aD"4^^(l) 
iaopl'^-^lf^r-l^l)  -oDM^^jd) 


m  =  0,1,2, 


From  the  definitions  of  f 


^--l)(t)  andfj-)(t)  it 


and  fi'^l(t)    it  follows  that  aD'"'"^'f^"':^^l)  = 

2  J  J- 


=  (7[a^J  and  that  aD™f  )"^.^  =  0  [a^"*"^]  and  consequently  that  the  denominator 
of  equation  (3.46)  does  not  vanish  when  a  is  sufficiently  small.   This  fact 
also  justifies  expanding  the  reciprocal  of  the  denominator  by  the  binomial 
theorem  in  powers  of  DV^"^^(l)/(iaaJD'™"^  If^'^^^^^l) ) .   In  this  way,  we 
find  that  for  m  >  1 


(3.U7)   c  =  I  h^°)(m)  + 


a 


2i(m-i) 


,    .  h^^^m)    i^h^°)(m)' 

Mm-M)h[°)(m)  ,^-^.-^^^^ 


haPb. 


Im 


3« 


[8h|°)(l)  +2h^°^(l)  +h^2)^^)j 


^  aca^j 


while  for  m  =  0  we  find  that 

(3.^)        c^  =  llia^h[°^0)  +  0[o^] 


,(o) 


,(2) 


where  the  quantities  h^   (m) ,  h^  ■^  (m)  are  defined  in  equation  (2.19) 
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Similarly,  we  find  in  the   case  of  H-mode   excitation  that  for  m  >  1 


(3.^9) 


c      =  2a^h^^^(m) 
m  o     ^ 


2ar 


'Ijn 


LiM^-l         3ni 


(^aS 


while  for  m  =  0 


(3-50) 


c     = 


h^^^O)   +|a^h|l)(0)    +6r[a3] 
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h.      Explicit  Results  for  Small  Values  of  a  =  ka 

U.l  General  Excitation.   In  what  follows  the  symbols  h^  ,  h^   ,  and 

(2) 
h^   denote  the  constants  appearing  in  equation  (2.I9).   We  distinguish 

between  E-mode  and  corresponding  H-mode  quantities  by  placing  an  overbar  on 

the  latter.   We  define  the  symbol  e  by 


(^.1)  e  =  \M   -  ^ 


The  transmission  coefficient  t  corresponding  to  an  incident  field 

m 

(E^'^^  H^"^^)  [see  equations  (2.U)  and  (2.5)]  maybe  defined  as  the  average 
power  radiated  to  infinity  in  a  period  on  the  shadow  side  of  the  screen 
divided  by  the  average  power  in  the  incident  wave  which  enters  the  shadow 
side  through  the  aperture  in  a  period.   It  is  easily  shown,  with  the  aid  of 
equations  (2.16)  and  (2.17),  that  for  m  >  1 

na!^  r"sino[|Af-1^0)|2  +  ]  cos  0  A^^^'^^  (O) -sin  0  A^  ^  (9)  |  ^ 


dO 


(^•2)    t   =       '° 


a  Re  /    /   E^  vv  H^   -e   rdrdo 
/   ~o    ~o   ~z      ' 


a|'^^(G)  and  A^'^"^^(0)  being  defined  as  in  equations  (3.IU)  and  (3.15)-  In 

equation  (^.2)   the  asterisk  denotes  complex  conjugation  and  e  denotes  the 

unit  vector  in  the  z-direction.   This  formula  applies  to  H-mode  as  well  as 
to  E-mode  excitation.   When  m  =  0  we  find  for  E-mode  excitation  that 

2nct    Tsin  0  |  cos^(j)^°  V^"^^0)cos  9  -  AJ°^9)sin  0|^+sin^<t>o°  ^2"^^  1^  }  ^^ 


(^.3)   t 


o  p2jt  pa 

a^- 


Re  f         P  E^°^/x  H(°).e   rdrdi 


S°)    ..  ^.^n=.^^  >.,.  l(°) 


The  same  formula  holds  for  H-mode  excitation  when  (|)^   is  replaced  by  (|)^ 
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In  this  subsection,  we  list  the  first  terms  in  the  power  series 

expressions  for  various  quantities  of  interest  -  general  E-  and  H-mode 

excitation  is  assumed.   (see  the  neighborhood  of  equation  (2.19)).   In  the 

next  subsection  we  specialize  these  general  formulas  to  the  case  of  plane 

wave  excitation  in  which  case  the  constants  h   ,  and  h  '^  ,  |-l  =  0,1,2,... 

j   =   0,1,2  are  given  by  the  expressions  listed  in  equation  (2.2^).   The  unit 

vectors  e  ,  ei ,  e^,  e  appearing  in  the  expressions  below  are  defined  at 

the  beginning  of  Section  2;   e  ,  e  and  e  denote  the  unit  vectors  parallel 
^  ^  '  ~x'  ~y     ~z 


to  the  X,  y  and  z  axes  respectively. 
(E-mode)   m  >  1 

In  the  aperture  CC  ,   we  find  that 


(h.k)    v[™^ 


rJ2) 


(o) 


sf:   i[(2m):J^(2m+l)[h^^^(m)+2h^°^m)Jr'"  e  a 

^^+1  m+1  ,  „,   ,  Is  „2/   ,  3n 
2    a   ml  r(m  +  -)  T  (m  +  ^) 


m+1 


r  m+Bl 


->-(T[a"^J, 


(^.5)    V, 


m-1 


m-1  m+1 
r   a 


2i(m+l)h[°''(m)  [m  e^  +  |  r^} 


m+1  ,  1. 

a    (m  -  g) 


+ 


mifh^^^(m)  +  2h^°''(m)] 
■  o  ^    ' o     -^ 

_  m-1,  2   Ix 
2a   (m  -  J-) 


-^,„0[a3]. 


where 
ih.6) 


Ira 


J- 


when  m  =  0 
when  m  f  0 


The  tangential  component  of  the  electric  field  in  the  aperture  is 
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.,    „.      ^        (-lfm(m+l)(2m-l)[(2m-2);]^r"^-^   e   g"^""^  ^ 

^      ''      ~t        ^Um-3     m+2.      ,v,     r-^,  ^    , 

2        -^  a        (m-l);  \/n  r(   ra  +  i  ) 


[h 


^^^(m)    +   (8ra+l|)h[°'' (m)    +  2h^°'' (m)  ]  f  sin  mcp^^   +  cos  mq?  e   } 


V2)^)    +2h^°)(m)\  ^2 


^°i  r  in 


+hi 


The  far  field  behavior   of  E  is    Fwe  remind  the  reader  that  H  =  e_  -a  E 
where  R  -^m.      See   equation   (2.l6)]: 


(U.8)      E 


(-l)"^./^i"^(sin  Q^-^^-"^ 
2^-^B.  r(m  +  |)r(m  +  |) 


"(m+l)(2m):h[°)(m) 


r(m  +  |) 


m 


(2m-2):fh^^^m)    +  2h^°^(m) 


+ 


r(m  +  |) 


I  sin  mcp  e_   +  cos   0  cos  mcp  e  1 


ikR 


m-cp-'      R 


+  ^.^r<x^j 


Im 


(R  ->«)) 


Finally,  the  transmission  coefficient  is 


-rt^(m+l)(2m+l)(2m-l):  [(2m-2):]^      /  2m+2 

^^'^^   *m  "  .  ,8m-liK  (o),  .|2-.1|,  ^   l.^li,   _^  3^  \M 
im2     |h^   (m)  I  P  (m  +  -)T    (m  +  ^)   \  (3 


[2m(2m-l)(m-fl)h[°)(m)-to(m  -  i)  (h^2)^2h^o)  ^^^  ^-12  ^  ^^'^Sm+Uj 
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When  m  =  0,   the  corresponding  results  are: 
(4.10)  v|°)    .  fg    [ka^l^[°ho)    -  h^2)(o)]    +  da?], 


f    .^  2iarh,^°^(0)  _ 

(u.ii)        v^-1)  =- — i — .a[cc3]. 


o^(o),  ,  (o)        2 

+8h^   ^(0)  sin  cp^       ^a  e 

(1..12)  E^  =  — ^— ^ 2 

rt  a   e 


5(a2-l)h[°)(0)-  2h^2)(o)j 


+  ^ cos  9^       ra  e^ 


+ 


OU"], 


(o)   ^  2^ 

n  a      e 


(4.13)  E  ~     ^^  /^"  ^      [cos   0[Ml-a2)J°^O)    +h^2)^Q)      ^^^      (o)^ 


-  4h$°)(0)    sin  cpi°)   .  ]    V     ^^[-^J  (R  -")' 


1      "-    '  ^o       ^^        R 


(1+.14)        t      = 


[|Ml-a2)h[°)(0)+h^2)^Q)|2^^^2^(o)^Q|^(o)^Q)|2^.^2^(o)^       ^2 
°  4^  |h$°)(0)|2   sin%(°°  ^"l^ 


l6a 


+ 


1 

r  4-1 


OM 


■3h- 


(H-Mode)   m  >  1 

In  the  aperture  (J^    we  find  that 

,  ,    Jn    r(2m):]2(2m+l)h^^^m)r™  €  a"'"'^ 

^^•^^^  ^1   =>-fl  ,  mfV,    1,^2,   ,i, +(^ta   J, 

2    m:a  T'Cm  +  -)r  (m  +  ^) 


r(l)/  \   m-l  m+1 
,   ^,    mh^  Mm)  r    a  ^     o 

2a    (m   -  1^) 
The  tangential  component  of  the  electric  field  in  the  aperture  is 

-im[(2m-2):j2(2rn-l)(m+l)  Vn  h^^^(m)    r"""^    e  a""""" 

'"^  ~t   "        o^m-3  m+2,      7~w  l^w      T~lTw      Tl^  ^ 

2        -'a        (m-l)ir(m   -  -)r(m  +  ^jrim  +  |) 

cos  mcp     e      +  sin  mcp        1  +  ^  )    e        +  (Tfa       1  , 


The  far  field  behavior  of  E  is 


,    ^ ^m    /-   .m-1    ,_      ^v    -(1),    N,    .      ^sm-1  2m+2 
_        (-1)     sjv:  1       m(2ra-2)lh^      (m)(sin  0)        a 

(^.18)        E ■ 

2^-^   ar(m   -  |)r(m  +  |)r(m  +  |) 


ikR 


[sin  mcpj^  eg  +  cos   9  cos  mcp^  |<p]V~  ^^  ^"^ 


The  transmission  coefficient  is 

jT^m(m+l)(2m+l)(2m-l):  [(2m)l]^|h^-^^(m)  1^  2m+2 

(^.19)    t   =  r: T-^ ^ •  "-  ,    s 

"^    i28-lh(°)(m)|2r2(m.|)r6(m.|)  (P^^^^ ) 


■X- 
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When  m  =  0  the   corresponding  results   are: 
(U.20)  4°)    =— ^^ -^    ff[a^]. 


(1..22)  E^  =  — ^- f— ^2 ^  -^    0[a5]  , 

3«a 


-8hS"^^(0)    cos^  cp^°''cos  Q   sin  0  a       ikE  ^^    „ 

(U.23)  E~— i ^^2 --—     +     Cr[aT  (R 

45rtia 


256aVi^^0)|^cc^ 
(ij-.2)      Plane  Wave  Excitation.      ¥e  wish  to  check  o\xr  results  for  m  =  0 


(^•2^)     t   =  ^ ^  ^  .  s ^   .  V  ^    +  (7(a^). 

o   in  TOR  2. |-(o),„^,2,-(o).*      ^^ 


and  m  =  1  against  those  collected  by  C.  J.  Bouwkamp  in  references  [3a]  and 
[313].   To  this  end,  we  note  that  Bouwkamps'  angle  G  corresponds  to  our 
angle  7  [see  the  neighborhood  of  equation  (2.21)]  and  that  Bouwkamps' 
plane  wave  is  the  linear  combination 

-k"^  cos  <j)  (E  ,  H  )  -  k"^sin  (t  (E  ,  H  ) 
"o  ~o  ~o  ^o^~o'  ~o 

of  the  plane  wave  fields  (E  ,  H  )  and  (E  ,  H  )  defined  in  equations  (2.21). 
■^  \~o'  ~o      ^~o'  ~o 

Setting  m  =  0  and  m  =  1  in  the  above  E-mode  and  H-mode  formulas,  inserting 
for  the  constants  h^"^'^(m)  and  h^"^'^(m)    (m  =  1,2)  appearing  there   the 
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expressions  in  the  right  members  of  equations  (2.25)  and 

finally  forming  the  appropriate  linear  combination  of  the  foior  fields 

thus  obtained,  we  arrive  at  the  results  listed  below. 

The  tangential  component  of  the  electric  field  in  the  aperture  is 


-2   cos   9      sin  Ore  q.  ., 

(^.25)  E,   =  ^°  °        -^        Q"^  ^ 


"t  T-r r—  a 


/2  2  -"       ,/2 

\/a     -  r  VS'  -^ 


2        12 
sin  i     cos  0   (a      -  -r   ) sin  c 
^o  o'  2 


i      f/    2        1  2.  ,,    2        1  2.     .    2^ 
+  cos  (p      '(a      -  ^r    )+(r      -  -a   jsm  0  cos 
"o   ^  2  2  o 


e 
~r 


/2     2 
8ia^ya   -r  12 

f-sin  (t  cos   9  cos  cp   +  cos  A   (1   -  -sin  9   )sin  cple 

a  L.     Tq     o    ^       'o^    2     o     ^-''=<p 


-^a[cc2J 


The  far -field  behavior  of  the  electric  field  is 

k^  3 
E  ~     -— ^  e„  /\  [(-Usin  i>  cos   9     +2   cos   A     sin  9     sin  9   sin  cp)e 
~  3jt     ~R        ^^  ^o  o  ^o  o  ^  ~x 

ikR 

_    Q 

+   (Ucos  (t   -  2  cos  (t  sin  9  sin  9  cos  cp)e  I    ~ 
~o  ~o      o  ~y-'    H 

The  fonnula  for  transmission  cross-section  is 

?  1; 
„   6ka  a      r  •  2i    2^  ^    2,1   2 ,   .  2_  -.  _^  ^r  6n 

where  the  transmission  cross-section  is  the  numerator  of  the  transmission 
coefficient. 

The  results  above  are  precisely  those  given  by  Bouwkamp  in  [3a, bj. 
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Appendix  I 
In  the  aperture  the  r  and  ^   components  of  the  scattered  magnetic 


field  H  are 

(m), 
H     =  cos  mi 


r  'm 


"P-  ^^  ^2  ^        ar  "^  Sr  \^    3z    JJ    ' 


(I.l) 

w               •         A  r  n2,jn-l  ^m„(m-l)  ,   m   Sv^™^ 

Hi  =   -sm  mo  (  kv„        +  -  V^        '    + =r- 

9                     '^m  I           2            r  r      o! 


r      dz       (  ' 


where 

(I  2)      r*^"^""""^    =  i  ^  frV^™'"'"^!    -  -  V^™'-^^ 
^  r   Str  \     2  /        r     2 

These  formulas  follow  directly  from  equations  (2.6)  and  (2-9) •   We  consider 

E-mode  excitation  and  H-mode  excitation  seperately 

E-mode  Excitation:   In  this  case  the  condition  that  H   =  H 
■  ~t   '-'o ,  t 

in  (JL   and,  equations  (I.l),  (I.2)  and  (2.4b)  imply 
^^•3^     r   o   -  ^  ^2    ^    ar    ^  Sr    az   ^ 


(1.1+)  -ik  -^  =  k^  v^-1^  -f  2i  r^"^-i)  +  S  _i 

or        2       r         r   oz 


Subtracting  equation   {l.h)   from   (1.3)   "we  get 


1_     ^  p(m-l)    _    .j^^(m) 
Laz 


=  0, 


from  which  we  infer  immediately  that 
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(1.8)        — ^ —  +  r^   '    -   ikV^  '    =   c  k(kr)  . 
^        '  dz  o      m 

where  c  is  an  arbitrary  constant.   Substituting  the  expression  of  Sv^   /5z  +  T^    ~   ' 
obtained  from  this  equation  into  equation  (I.U)  and  combining  terms  in  an 
obvious  manner  we  find  that 

(1.7)         yX  =  ~T  '^        ^r-  r  v^     -  mc  (kr) 

2        k     or  V    o  /     m 

Substituting  in  turn  this  expression  for  V^     into  the  right  member 
of  equation  (1.2),  putting  the  resulting  expression  for  P      into  equation 


(1.6)  and  noting  that 


1  a  /    K 


(1-8)  Zh[r 


(m)  \         _  o,  a2v(™) 


o 


k 


m 


2- 


az 


r   Sr  \        Bz 
[see  equqtion   (2.Ud)J   we  get 

oz  dz 

Equations  (l.?)  and  (1.9)  will  be  recognized  as  the  conditions  listed  in 
equation  (2.1^+). 

H-mode  Excitation:   Employing  equations  (l.l),  (1.2)  and  (2.5b), 
we  obtain  the  analogs  of  (I.3)  and  (1.4)  [See  footnote  7  J 


rr  in^       ^  (   o   1    v2-(m-l)  ^  S  {   ^\ 


(111)    ^--k^iv(-^)   5  iX^ 

^^•^^^      Sz   -   ^  m  ^2     -  Sz   -    ¥ 
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where 

~rarV2         J    '  r     2 

Differentiating  equation  (l.ll)  with  respect  to  r  and  adding  the  resulting 
equation  to  (I.IO)  we  get 

(1.12)  |:(-^^""'^4""'^  =° 

This  equation  implies 

(1.13)  V^'-^=  mc  (kr)"^-^ 

where  c  is  an  arbitrary  constant.   Substituting  this  expression  for 

m 

-(m-1)  ^^^^  equation  (l.ll)  we  find  that 

il.lk)  V  "^^=  -  ^  -  c  k(kr)"^, 

^    '         1      dz      m      ' 


Equations  (1.13)  and  (I.l4)  will  be  recognized  as  those  listed  in 
(2.15). 
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FOOTNOTES 

1.  Marcuvitz's  solution  is  outlined  in  our  appendix  of  Bouvkamp's 
report. 

2.  K.  A.  Lur'e's  paper  is  an  outgrowth  of  a  work  by  N.  W.  Lebedev  and 
I.  P.  Skal'skaya  [lo]  who  considered  the  special  case  of  normal 
incidence  of  a  plane  wave  on  a  disc.   Both  Lur'e  and  Williams  are 
concerned  with  the  diffraction  by  a  disc.   The  solution  of  the  disc 
and  aperture  problems  are  of  coiurse  related  through  Babinet's 
principle. 

3.  We  are  grateful  to  W.  D.  Collins  for  bringing  to  our  attention  the 
existence  of  yet  another  solution  of  the  vector  aperture  problem  by 
J.  P.  Benkard  University  of  Michigan  thesis  entitled  "On  an  axi- 
symmetric  vector  boundary  value  problem  with  specal  reference  to  the 
circular  disc." 

k.        Representations  of  the  kind  appearing  in  equation  (1.1)  were  originally 
introduced  by  J.  Bazer  and  A.  Brown  [lo]  for  the  case  m  =  0,  by  W.  D 
Collins  [llj  for  the  cases  m  =  0,  1  and  by  J,  Bazer  and  H.  Hochstadt 
[6j  for  general  values  of  m  >  0. 

5-    These  integro-differential  equations  lead  to  Fredhom  integral  equations  of 

the  second  kind  having  Jones-type  kernels  [see  reference  [l3a]] .   Such 

integral  equations  are  useful  in  establishing  the  existence  of  f^^'^'^(t) 
/   -i  \  _  1 

and  f  (t);  however,  as  explained  in  reference  [j] ,  they  do  not 
provide  an  effective  method  for  calculating  f-j"^  and  fp  ~  when  m 
is  large  even  when  ka  is  small. 

6.  It  would  be  tedious  but  not  impractical  to  obtain  a  few  additional 
terms. 

7.  Wherever  E  and  H-modes  are  discussed  together  and  there  is  a  risk  of 
confusing  E-mode  quantities  with  H-mode  quantities,  we  employ  an  over- 
bar  to  distinguish  the  latter  from  the  former  as  in  equations  (2.4), 
(2.5)  and  (2.l4)-(2.15). 
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8.  Various  formulations  of  the  edge  conditions  for  the  vector  problem 
will  he  found  in  the  papers  by  C.  J.  Bouwkamp  [^Sa^ti"!  J.  Meixner  pLbl 
and  D.  S.  Jones  [iStiJ  .   We  employ  here  a  formula  due  to  N.  Marcuvitz 
I^See  appendix  of  reference  [_3t)]l  • 

9.  This  assumption  is  employed  to  merely  simplify  the  general  fonnulas 
for  the  transmission  coefficients  j^see  equations  {k.9    ),(h.lk)    and 
(U. 19) , (4.24)J  ;   it  is  not  used  in  any  other  calculation. 

10.  This  section  is  patterned  after  Sections 2  and  3  of  J.  Bazer  and 
H.  Hochstadt  [7] . 

11.  'See  property  (i)  below. 

12.  When  j  =  2  and  m  =  -1,  we  require  that  (ii)  be  satisfied  for  (ka)fp   (t). 

13.  Property  (2)  (f^"^^(t)   is  regular  in  the  domain  |t|  <  1  +  A/a,  A  >  0) 
is  used  in  this  connection. 

ik.        The  limit  of  the  integral  term  as  z  approaches  zero  is  zero  owing  to 
the  oddness  of  f|"^(t). 

15.  Equations  (3.2O) -(3.2i+)  hold  for  all  m  >  0  except  that  the  factor 
l/a  in  the  right  members  of  (3.22)  and  (3-24)  must  be  replaced  by 
l/a  when  m  =  0. 

16.  Hereafter,  we  shall  use  the  abbreviation  a  =  ka. 

17-   We  ignore  for  the  moment  the  presence  of  the  undetermined  constant  c  . 

m  h 

18.  When  m  =  0  the  factor  a  in  the  right  member  must  be  replaced  by  a  . 

19.  This  restriction  on  the  degree  of  B   (t)  may  be  verified  by  substituting 
f   (n/t])  into  equation  (3.  33)  ,  expanding  in  power  of  a  and  comparing 
coefficients  of  like  powers. 

20.  The  radius  of  convergence  of  this  power  series  has  not  been  determined. 
For  the  case  of  plane  wave  excitation  and  m  =  0  it  is  known  to  be  at 
least  1/2  [see  reference  [15J] • 

21.  6.  .  denotes  the  Kronceker  delta:   8.  .  =  0  if  r  j^  j  and  S.  .  =  1  if 

i  =  J- 

22.  ¥e  have  multiplied  the  numerator  and  denominator  of  this  fraction 
by  a  =  ka  to  avoid  having  to  discuss  the  case  m  =  0  separately. 
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